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Electromagnetic  Penetration  into  a  Finite  Coaxial 


Cylinder  with  a  Recessed  Inner  Conductor 


Larry  Rispin  and  David  C.  Chang 

\ 

Abstract 

-4 

The  electromagnetic  penetration  into  a  finite  coaxial  cylinder, 
illuminated  by  a  uniform  plane  wave  of  arbitrary  incidence  is  investi¬ 
gated  by  combining  the  results  of  several  individual  Wiener- Ho  pf 
analyses.  The  circular  to  coaxial  waveguide  junction  within  the  cylinder 
is  fully  characterized  through  such  a  technique.  A  simple  equivalent 
circuit  far  this  junction  and  the  other  elements  of  the  system  are 
combined  into  an  equivalent  Norton  source  and  equivalent  load  admittance 
which  is  valid  under  most  practical  situations.  The  power  dissipated 
within  the  coaxial  line  section  of  the  system  is  then  easily  calculated 
far  any  given  load,  except  those  in  the  vicinity  of,  the  "perfectly 
matched"  load. 
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1.  Introduction 


The  study  of  the  electromagnetic  penetration  into  cylindrical 
enclosures  is  of  great  practical  concern  due  to  the  effects  it  may  have 
on  the  mechanisms  or  circuitry  contained  within  the  enclosure.  The 
most  prevalent  case  of  this  high  frequency  electromagnetic  penetration 
takes  place  through  intentional  as  well  as  unintentional  apertures  on 
the  enclosure,  which  are  formed  by  cracks,  seams,  access  doors,  etc. 

The  study  of  this  type  of  aperture  coupling  is  well  documented,  [ 1 ] - 
[10].  Another  form  of  electromagnetic  coupling  into  a  cylindrical  en¬ 
closure  occurs  through  the  open  end  of  a  circular  cylinder.  This  has 
been  considered  in  the  particular  case  of  a  finite-length  thin  cylinder 
open  at  one  end  and  illuminated  by  a  uniform  plane  wave  by  Chang,  Lee 
and  Rispin  in  [11].  This  investigation  showed  the  electromagnetic 
penetration  into  the  open  end  of  the  cylinder  to  be  predominantly  in 
the  form  of  evanescent  circular  waveguide  TMQN  modes.  Since  these 
waveguide  modes  are  cut  off,  their  effects  at  a  sufficient  distance 
away  from  the  end  of  the  cylinder  are  usually  quite  small.  However, 
if  an  insulated  conductor  within  the  cylinder  is  near  the  open  end, 
the  evanescent  waveguide  modes  may  couple  into  a  TEM  mode  supported  by 
the  insulated  conductor  and  the  cylinder.  Energy  may  then  travel 
practically  unattenuated  along  this  two-conductor  system  to  delicate 
circuitry  within  the  structure.  In  order  to  develop  a  quantitative 
feeling  for  this  type  of  electromagnetic  coupling,  the  theoretical  model 
shown  in  Figure  1  was  proposed.  The  two-conductor  transmission  line  is 


simulated  by  a  coaxial  line  whose  center  conductor  is  recessed  from 
the  end  of  the  cylinder.  A  load  admittance  is  located  at  the  opposite 
end  of  the  center  conductor  which  simulates  the  lumped  admittance  of 
the  circuitry  in  the  practical  situation.  For  convenience  we  assume 
the  cylindrical  enclosure  to  be  open  at  only  one  end  and  capped  by  a 
metal  plate  at  the  other. 

The  analysis  of  the  model  shoim  in  Figure  1  may  be  conveniently 
divided  into  three  major  parts.  First  the  initial  penetrating  currents 
at  the  end  of  the  cylinder  may  be  determined  through  the  Wiener-Hopf 
technique  as  demonstrated  by  Chang,  Lee  and  Rispin  [11].  The  circular 
waveguide-to-coaxial  waveguide  transition  must  be  characterized  next. 
This  characterization  forms  the  bulk  of  this  report,  namely  Sections  2, 
3,  and  4.  And  finally  simple  transmission  line  theory  may  be  used  for 
the  final  section  between  the  ci rcular-to-coaxial  junction  and  the 
load  admittance.  In  this  manner  we  may  determine  the  load  current  and 
power  due  to  the  incidence  of  a  uniform  plane  wave  of  known  intensity 
and  angle  of  incidence. 

As  mentioned  above,  the  major  part  of  this  report  is  devoted  to 
the  complete  characterization  of  the  circular-to-coaxial  waveguide 
transition.  We  note  that  in  the  case  of  an  incident  wave  propagating 
toward  the  transition  from  the  coaxial  region,  (instead  of  incident 
waves  propagating  toward  the  transition  in  the  circular  waveguide  region 
which  are  due  to  the  plane  wave  impinging  upon  the  cylinder),  the  same 
structure  also  characterizes  the  problem  of  a  typical  coaxial  open 
circuit  termination,  [12]  and  [13].  Thus  in  Section  2  a  "spectral" 
equation  and  an  auxiliary  equation  in  the  Fourier  transform  domain  are 


i 


developed  for  the  transition  region.  Section  3  defines  the  solution 
of  the  circular-to-coaxial  waveguide  problem  in  terms  of  a  "composite" 
problem  with  waves  incident  from  both  sides  of  the  transition.  In 
Section  4,  the  Wiener-Hopf  technique  is  utilized  to  obtain  the  scattered 
currents  from  the  transition.  Equivalent  circuit  representations  of 
the  circular-to-coaxial  waveguide  junction  are  given  in  Section  5.  And 
finally  in  Section  6,  the  characterization  of  the  circular-to-coaxial 
waveguide  transition  of  the  previous  sections  is  applied  to  the  system 
shown  in  Figure  1. 
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2.  Formulation  of  the  Spectral  liquation 


Consider  the  coaxial  system  shown  in  Figure  2,  in  which  we  have 
a  perfectly  conducting  surface  at  p  =  b  concentric  about  an  arbitrary 
(perfectly  conducting  or  nothing)  cylindrical  surface  of  radius,  a. 
Leaving  the  surface  at  p  =  a  arbitrary,  allows  the  latter  specification 
of  it  as  being  either  an  infinitely  long  perfectly  conducting  cylinder, 
or  a  semi-infinitely  long  perfectly  conducting  cylinder,  or  no  conducting 
surface  at  all.  We  assume  an  exp(-io)t)  time  variation  and  define  the 
Fourier  Transform  pair  as. 


F(cO 


■  f 


Cl  ^  10tz  . 

F(z)e  dz 


(2.1) 


F(z) 


2-n 


(2.2) 


Considering  only  axially  symmetric  fields  and  currents  allows  us 
to  write  the  transform  of  the  z-directed  electric  field  as, 

{C2AJQ(^P)  ;  0  <  p  <  a 

(2.3) 

C2[BJ0Up)  +  CY0(i;p)]  ;  a  <  p  <  b 

where  the  unknown  constants.  A,  B,  and  C,  are  determined  by  imposing 
the  appropriate  boundary  conditions,  and, 

C  =  /k2  -  a2  =  i/ a2  -  k2  ;  i.e.,  I  (£)  >  0 

’  mv  '  - 

for  all  a  (2.4) 


The  other  field  quantities  under  the  assumption  of  axial  symmetry 
may  be  given  in  terms  of  Ez(p,a)  as, 


Figure  2.  A  coaxial  system  consisting  of  a  perfectly  conducting  outer 
conductor  of  radius  b  with  an  arbitrary  (perfectly  con¬ 
ducting  or  no  conductor  at  all)  coaxial  inner  surface  of 
radius  a  . 
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Ep(P.“)  ■  -  77  -JJ  EjCP.O)  np(p,a)  =  0  (2.5) 

^(P.P)  -  0  n  CP.o.)  -  -£■  £  E  (p,a) 

?n 

(2.6) 

Hz(P,a)  (2.7) 

For  the  system  described  in  Figure  2,  the  ap  ; •H ate  transformed 

boundary  conditions  are. 


Kz(b,a) 

=  0 

(2.8) 

„  + 

Ez(a  ,a) 

=  Ez(a  ,a) 

(2.9) 

~  + 

}(<*) 

H0(a  ,o) 

-  HA(a  ,a)  -  “ 

4>  2-rra 

(2.10) 

+ 

-  . 

H^(b  ,oi) 

-  Vb  ■“>  “  2,b 

(2.11) 

where  (2. 8)specifies  that  the  z-directed  electric  field  at  p  =  b  must 
be  zero,  (2.9)  enforces  the  continuity  of  the  z-directed  electric  field 
at  p  =  a,  (2.10)  expresses  the  discontinuity  in  the  ^-directed  magnetic 
field  at  p  =  a  in  terms  of  a  z-directed  current  there,  and  (2.11) 
expresses  the  discontinuity  in  the  ^-directed  magnetic  field  at  p  =  b, 
(note,  H(h+,o)  =  0)  in  terms  of  a  z-directed  current  there. 

The  constants  of  (2.3)  may  be  determined  in  terms  of  I  (a)  by 

a 

utilizing  the  boundary  conditions  in  (2.8),  (2.9),  and  (2.10).  With 
the  constants  of  (2.3)  in  this  form,  wo  may  take  p  =  a  yielding  the 


"spectral"  equation, 


Ez(a,o)  =  ~  g  ;2M(a)N(a)Ia(a)  (2.12) 

where, 

M(a)  =  ir[J0(ea)Y0tth)  -  J0(4b)  YQ(Ca)  ]  (2.13) 

(Ca) 

Nlo)-3^KT  (2-1J) 

(2.12)  is  the  equation  utilized  in  the  Wiener-Hopf  analysis  for  the 
circular  waveguide-to-coaxial  waveguide  transition  to  be  discussed  in 
Section  4. 

Alternatively  the  constants  of  (2.3)  may  be  determined  in  terms  of 
Ib(a)  by  utilizing  the  boundary  conditions  in  (2.8),  (2.9),  and  (2.11). 

With  this  new  set  of  constants  and  for  p  =  a  in  (2.3)  we  find  the  relation¬ 
ship  between  the  transforms  of  the  currents  on  the  p  =  a  and  b  conductors, 

Ib(oO  =  -N(a)la(a)  (2.15) 

(2.15)  will  alow  us  to  eventually  solve  for  the  scattered  current  on 
the  p  =  b  conductor  for  all  z  from  only  one  equation. 
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3.  Definition  of  a  Composite  Problem  for  the  Circular  Waveguide- to- 
Coaxial  Waveguide  Transition 

3.1  The  individual  systems 

To  characterize  the  circular  waveguide-to-coaxial  waveguide  transi¬ 
tion  we  consider  the  three  systems  shown  in  Figure  3.  Figure  3B  illus¬ 
trates  an  infinitely  long  circular  waveguide  of  radius,  b,  with  an 
arbitrary  source  at  z  =  -L.  The  incident  current  on  the  p  =  b  conductor 
for  z  >  -L  is  given  by, 


.  B 


(2)  =  I 

N=1 


(Ib^ON  C 


-v 


•L  <  z  < 


(3.1) 


where  (Ib)0N  i-s  the  coefficient  of  the  TMQN  circular  waveguide  mode 
current  (at  z  =  0)  and  iybN  is  the  propagation  constant  of  the  ™0N  mode 
given  by, 


(3.2) 


and  pN  is  the  Nth  ordered  zero  of  Jq(x),  i.e., 

J0(pN)  =  0  ;  N  =  1,2,3,  ...  (3.3) 

Associated  with  this  incident  current  is  an  incident  z-directed  electric 
field  at  p  =  a  given  by, 


rB  .  r*  bN 

-  [  (£,)„„  c  ;  -L  <  a  <  - 

N=1 


(3.4) 


B  B 

The  model  coefficients  of  Ez(a,z)  and  I  (z)  are  easily  shown  to  be 


related  by, 


Figure  3.  (A)  Semi-infinite  circular  waveguide  to  semi-infinite 
coaxial  waveguide,  sources  at  z=-L  and  +L. 

(B)  Infinite  circular  waveguide,  source  at  z  =  -L. 

(C)  Infinite  coaxial  waveguide,  source  at  z  =  +L. 
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^  Von 


(Von 


(3.5) 


Also  it  is  obvious  that  no  current  exists  at  p  =  a,  thus, 

lJU)  =0  ;  -L  <  z  <  °"  (3.6) 

In  Figure  3C  an  infinite  coaxial  waveguide  with  an  inner  conductor 
of  radius,  a,  and  outer  conductor  of  radius,  b,  is  shown  with  an  arbitrary 
source  located  at  z  =  L.  The  incident  current  for  z  <  L  is  given  by 


■b<z>  ■  I  (>h> 

M=0 


C) 

bJ0M 


W 


-<a  <  2  <  L 


(3.7) 


where  (I^)QM  is  the  coefficient  of  the  TMQM  coaxial  waveguide  mode 
current  (at  z  =  0)  and  iYcM  is  the  propagation  constant  of  the  TMQM  mode 
given  by, 


(3.8) 


And  p^,  which  is  a  function  of  the  ratio  a/b,  is  the  Mth  ordered  zero 
of  the  cross  product  of  JQ  and  Yq,  i.e.. 


J0(vpcM)Y0(fW 


VpcM^VvpcM^  =  0  »  M  =  1’2,3’ 


(3.9) 


where 


v  =  a/b 


(3.10) 


And  the  special  case. 


Pc0  =  0  (3.11) 

corresponds  to  the  TMQ0  or  TEM  mode. 

Associated  with  the  current  on  the  p  =  b  conductor  is  a  current 
on  the  p  =  a  conductor  given  by. 


oo 


1 

M=0 


t-VoM 


<  z  <  L 


(3.12) 


The  modal  coefficients  of  the  currents  on  the  inner  and  outer  conductors 
are  related  by, 

<‘S>0N  *  '3-t5> 

And  finally  for  the  system  of  3C  it  is  obvious  that. 


Ez(a,z)  =•  0  ;  -<*>  <  z  <  L  (3.14) 

The  system  to  be  investigated  is  the  circular  waveguide-to-coaxial 
waveguide  transition  shown  in  Figure  3A.  The  region  for  z  <  0  is 
identical  (including  the  source  at  z  =  -L)  to  that  of  the  circular  wave¬ 
guide  in  Figure  3B.  And  the  region  for  z  >  0  is  identical  (including 
the  source  at  z  =  L)  to  the  coaxial  waveguide  of  Figure  3C.  The  field 
quantities  in  Figure  3A  may  be  expressed  in  terms  of  the  incident 
quantities  of  Figure  3B  and  3C  plus  scattered  quantities  due  to  the 
transition  at  z  =  0,  i.e.. 


3 


E®(a,z)  +  E*  (a.z) 


Ez(a,z)  = 


l*(z)  - 


‘atz>  *  '.**<*> 


and, 


=  -( 


‘‘w  *  'b'W 


;  -L  <  z  <  0 
;  0  <  z  <  L 

;  -L  <  z  <  0 

;  0  <  z  <  L 

;  -L  <  z  <  0 

;  0  <  z  <  L 


(3.15) 


(3.16) 


(3.17) 


where  the  and  signs  on  the  superscript  "s"  designate  the  scattered 
term  for  the  regions,  z  >  0  and  z  <  0,  respectively. 

Since  the  inner  conductor  in  Figure  3A  would  also  carry  circular 
waveguide  TM^,  (*■  =  1,2,3,...),  mode  currents  on  its  internal  walls  in 
addition  to  the  coaxial  waveguide  TMq^  mode  currents  on  the  outer  wall, 
it  is  convenient  at  this  time  to  define  the  propagation  constant  for 
these  modes. 


where  has  the  same  definition  of  in  (3.3). 

An  additional  comment,  in  the  mode  current  definitions  of  this 


section  we  have  adopted  the  convention  of  indexing  incident  mode  currents 
and  fields  with  capital  letters,  M  and  N,  etc.,  and  scattered  mode 
currents  and  fields  with  lower  case  letters,  m,  n,  and  t,  etc..  This 
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convention  will  be  used  throughout  this  report  and  should  be  understood 
in  Appendix  B  concerning  the  discussion  of  so-called  "wel 1 -coupled"  modes. 

3.2  Construction  of  the  composite  problem. 

The  spectral  equation  in  (2.12)  and  the  supplementary  equation  in 
(2.15)  were  derived  under  the  assumption  that  all  of  the  implicated 
field  quantities  possessed  Fourier  Transforms.  Unfortunately  if  we 
move  the  sources  at  z  =  -L  and  z  =  +L  to  -°°  and  +»,  respectively, 

(i.e.,  L  -*■  «°),  the  incident  fields  and  currents  in  (3. 15)- (3. 17)  which 
do  not  satisfy  the  radiation  condition  of  an  outgoing  wave,  will  not  in 
general  have  Fourier  transforms  valid  in  the  "analytic  strip".  This 
analytic  strip  is  defined  as  -Im(k)  <  Im(a)  <  Im(k)  as  depicted  in 
Figure  4,  (note  as  usual  in  a  Wiener-llopf  analysis  we  assume  a  slightly 
lossy  medium,  i.e.,  k  =  k'  +  ik"),  and  contains  the  inverse  Fourier 
transform  contour  in  (2.2)  designated  as  rQ.  In  order  to  circumvent 
this  difficulty,  we  propose  the  simultaneous  solution  of  three  super¬ 
imposed  problems  which  features  the  standard  inverse  Fourier  transform 
contour,  Tq,  in  Figure  4.  By  combining  the  field  quantities  of  the 
three  systems  in  Figure  3  as, 

"Composite"  =  #3A  -  #3B  -  #3C  (3.19) 

the  composite  z-directed  electric  field  at  p  =  a  is  given  by, 

Ez(a,z)  =  E*(a,z)  -  E®(a,z)  -  E^(a,z)  (3.20) 

The  expressions  for  the  other  composite  quantities  are  similar  in  form 
to  (3.20).  Consequently  we  have, 
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Ez(a,z) 


Vz) 


Vz> 


f  Ez"(a,z) 


< 


-L  <  z  <  0 

0  <  z  <  L 


-L  <  z  <  0 

0  <  z  <  L 


-L  <  z  <  0 

0  <  z  <  L 


(3.21) 


(3.22) 


(3.23) 


The  composite  quantities  in  (3.21)  through  (3.23)  may  be  interpreted 
as  the  "extended"  scattered  (z-directed)  electric  field  at  p  =  a  and 
the  "extended"  induced  currents  on  the  p  ■  a  and  p  *  b  conductors, 
respectively,  for  the  system  in  Figure  3A.  The  term,  "extended",  refers 
to,  in  the  case  of  I  (z)  in  (3.22),  the  extension  of  the  induced  current 

a 

at  p  ■  a  from  the  coaxial  region  into  the  circular  waveguide  region. 

This  extended  induced  current  combined  with  the  incident  current  yields 
the  necessary  result,  that  the  total  current  at  p  =  a  for  z  <  0  is  zero. 
Ez(a,z)  in  (3.21)  and  1^ (z)  in  (3.23)  may  be  similarly  interpreted. 

These  "composite"  or  "extended"  quantities  satisfy  the  radiation  condi¬ 
tion  of  out-going  waves  and  may  be  used  in  the  spectral  equation  in 
(2.12)  and  the  auxiliary  equation  in  (2.15).  Removing  the  sources  of 


17 


the  composite  system  in  (3.19)  to  an  infinite  distance  away  from  the 
transition  at  z  *  0,  (i.e.,  L  -*■  °°),  does  not  inhibit  the  Fourier 
transformation  of  the  composite  terms  in  (3.21)  through  (3.23). 


4.  Wiener-Hopf  Solution 

Allowing  L  -*■  *  in  (3.21)  through  (3.23)  and  taking  the  Fourier 
transform  (2.1)  of  each  quantity  results  in  the  following  transformed 
composite  terms, 


Ma.a)  =  E*  (a, a)  -  £  (E?) 

N=1 


B) 

z)  ON 


(°+iYbN) 


(4.1) 


I  («)  -  S  +  («)  -  l  U?) 


M=0 


a  OM 


i 

(o‘iYcM) 


(4.2) 


Vo) 


I  CiKB) 


N=  1 


b  ON 


_ i 

(a+i 


iYbN} 


(4 


The  "+"  and  subscripts  appearing  in  (4.1)  and  (4.2)  signify  that  the 
function  is  analytic  in  either  the  upper  half  (-Im(k)<Im(a)<“>)  or  the  lower 
half  (-•<  Im(a)  <Im(k) ) ,  respectively,  of  the  transform  plane  shown  in  Fig.  4. 
We  note  in  passing  that,  in  order  to  include  incoming  waves  from  the 

region  P  <a  and  z  >  0  we  need  only  to  insert  the  term 


into  (4.2),  (1°)  .  being  the  coefficient  of  the 

3  U  L 

r  th 


mode  at  z  *  0.  The  regions  of  analyticity  for  the  terms,  (a, a) 
and  I  (a)  were  determined  using  the  procedure  described  by  Mittra  and  Lee 

ft* 

in  Section  3-3  of  [14]  with  the  known  asymptotic  behaviors. 
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and 
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as  z  -  00 


J.S+,  ,  ikz  .  m 

Ijj  (zj  e  ;  as  z  +  +  00 


(4.6) 


(4.7) 


the  same  procedure  reveals  I^fa)  to  be  analytic  within  the  analytic 
strip,  -Im(k)  <  Im(a)  <  In(k)  . 

Returning  to  the  spectral  equation  (2.12),  which  will  form  the  basis  of 
the  Wiencr-llopf  analysis,  we  express  the  factorization  of  the  kernels 
(see  Appendix  A)  as, 


M(a)  =  M+  (a)  M_(a) 

(4.8) 

and 

N  (a)  =  N+(a)  N_(a) 

(4.9) 

The  subscripts  "+"  and  indicate  the  region  of  analyticity  (upper 

half  plane  and  lower  half  plane,  respectively)  of  the  factor.  Upon  the 
substitution  of  (4.1),  (4.2),  (4.8),  and  (4.9)  into  (2.12)  and  after 
some  rearrangement  we  arrive  at  a  form  of  the  Wiener-Hopf  equation  given  by 


(k-o)M_(o)N  (a) 

-  ^  X  (k+a)M+(a)N+(a) 


(4.10) 


Following  a  typical  Wiener-Hopf  decomposing  procedure,  the  mixed  functions 
in  (4.10)  may  be  written  in  the  form. 


E  (a.al  °° 

i  *  '  r.  _ ^~z^0N _ 

(k-a)M_ (a)N_ (a)  (k+iYbN)M+(iYbN)N+(iy 


r  i  l 

bN1  L  J 


+; 


21 


(MK.WH.Wi.W 


CO 

'Ik  M=0(‘Ia-l0M(-k+1YcM^M+(-lYcM^N-*-(-:iTcM->  [  (a-iY„u)] 

r  (lfloN _  r  i  1 

+  N=1  tk+iYbN)M.(iYbN)N.(iYbN)  LC«-iYbN) J 
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(4.11) 


where  we  have  used  the  properties  from  Appendix  A, 

M_(-a)  =  M+ (a)  (4.12) 

N  (-a)  =  N+(a)  (4.13) 

and  the  definitions  of  the  composite  field  E  (a,a)  and  current  I  (a) 
as  given  in  (4.1)  and  (4.2),  respectively.  Now  since  the  LHS  of  (4.11) 
is  analytic  in  the  lower  half  plane  and  the  RHS  of  (4.11)  is  analytic  in 
the  upper  half  plane,  one  can  show  by  the  use  of  Liouville's  Theorem  [14] 
that  both  sides  are  equal  to  an  entire  function,  P(a),  which  is  analytic 
for  all  a  . 

From  the  edge  condition  it  is  easily  shown  that 

Is  (a)  ~  0(a"*)  as  jaj  -*■  00  in  the  upper  half  plane  (4.14) 

From  Appendix  A, 

M+(a)N+(o)  'v  0(a ~^)  as  |a|  -+■  «*  in  the  upper  half  plane  (4.15) 


Using  (4.14)  and  (4.15)  in  the  KHS  of  (4.11)  yields 

P(a)  *  0  for  all  a  (4.1b) 

Consequently  we  obtain  from  (4.11), 
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The  use  of  (2.15)  and  (4.18)  and  the  inverse  Fourier  transform  in 
(2.2),  now  provide  an  explicit  expression  for  the  composite  current  on 
the  p  =  b  conductor,  which  is  given  by 
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(4.19) 


where  the  contour  rQ  is  shown  in  Figure  4. 

Ifa(z)  for  z  <:  0  is  obtained  by  deforming  rQ  upward  enclosing  the 
poles  (from  the  integrands  of  (4.19))  at  a  =  iYcM  and  the  poles  of  N  (a) 
at  a  =  iy^.  The  residues  of  these  latter  poles  are  easily  determined 
by  making  the  substitution,  N  (a)  =  N(a)/N+(a)  from  (4.8),  the  end 
result  being 
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A  comparison  of  (3.23)  and  (4.20)  reveals  the  common  incident  term,  -I^(z), 


leaving  the  remaining  terms  in  (4.20)  to  be  identified  as  the  scattered 
currents;  the  sum  over  N  being  the  scattered  currents  due  to  the 


incident  circular  waveguide  TM  mode  currents  from  z  <0  and  the  sum 


over  M  is  the  scattered  currents  due  to  the  incident  coaxial  waveguide 


TM...  mode  currents  from  z  >  0. 
(JM 


To  find  the  composite  current,  I^fz),  for  z  zl  0  we  return  to  (4.19) 


and  deform  the  contour,  downward  enclosing  the  poles,  a  =  -k  and 


,-1 


-iy.M,  and  the  poles  of  [M  (a)]  located  at  a  =  -iy  -  The  residues  of 
dN  +  cm 

these  latter  poles  are  easily  determined  by  making  the  substitution, 

M  (a)  =  M(a)/M  (a),  the  end  result  being 
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A  comparison  of  (3.23)  and  (4.21!  immediately  reveals  the  common  incident 

g 

term,  -I^(z).  The  remaining  terms  in  (4.22)  are  thus  the  scattered 
currents  due  to  the  circular-to-coaxial  waveguide  transition  at  z  =  0. 

We  note  in  passing  that  N2 CiYcm)  =  Jo (VpCm^J0 ^cnr*  *  for  m  *  0  l15l> 
that  the  term  [1-  N(iYcm)]  1  is  always  finite  for  m  >_  1.  The  situation 
in  which  a  mode  in  system  B  of  Figure  3  (i.e.,  eigenvalues  of  ± iY^n) > 
approaches  the  same  distribution  of  a  mode  in  system  C  of  Figure  3, 

(i.e.,  eigenvalues  of  iiy^)  is  discussed  in  Appendix  B. 


so 
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5.  Equivalent  Circuits  for  the  Circular-to-Coaxial  Waveguide  Junction 

The  circular-to-coaxial  waveguide  junction  may  be  described  in 
terms  of  a  current  scattering  matrix  representation,  i.e.. 

[I']  =  [C]|I+]  (5.1) 

where  [I+]  is  the  incident  current  vector,  [I  ]  is  the  reflected 
current  vector,  and  [C]  is  the  current  scattering  matrix  which  contains 
the  interrelationships  between  the  various  incident  and  reflected  currents. 
Since  an  infinite  number  of  modes  may  exist  on  both  sides  of  the  junction 
the  terms  in  (5.1)  are  in  theory,  infinitely  dimensional.  However  in  most 
practical  situations  only  the  lowest  order  modes  [most  dominant  modes)  need 
to  be  considered. 

The  simplest  two-port  version  of  (5.1)  occurs  when  only  the  most 
dominant  modes  in  each  region,  the  circular  TM^  and  the  coaxial  TMQ0  (TEM) 
modes,  are  retained.  This  reduction  requires  a  demonstration  showing  that 
the  effects  of  the  higher  order  modes  are  indeed  negligible.  To  proceed, 
we  assume  that  these  higher  order  mode  effects  are  minimal  compared  to 
the  retained  dominant  modes,  and  offer  the  equivalent  circuit  shown  in 
Figure  5a.  For  this  equivalent  circuit,  (5.1)  is  reduced  to 
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in  which  we  have  arbitrarily  designated  the  circular  TM^  mode  currents 
(incident  and  reflected)  with  the  subscript  "1"  and  the  coaxial  TMqq(TEM) 
mode  currents  with  the  subscript  "2".  The  current  scattering  parameters  in 
(5.2)  are  explicitly  known  from  (4.20)  and  (4.21)  to  be: 
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Figure  5.  Equivalent  circuits  for  the  dominant  modes  in  the 
circular-to-coaxial  waveguide  junction. 


(a)  Two-port  current  scattering  parameter  representation. 

(b)  Two-port  admittance  parameter  representation. 

(c)  One-port  discontinuity  capacitance  representation. 
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In  obtaining  (5.3)  through  (5.6)  we  have  replaced  YcQ  by  -ik,  which 
is  a  result  from  (3.8)  and  (3.11).  The  current  scattering  parameters  in 
(5.3)  through  (5.6)  can  be  shown  to  be  those  of  a  lossless  reciprocal 
junction  as  would  be  expected  on  the  basis  of  physical  reasoning. 

A  more  conventional  equivalent  circuit  representation,  using  admittance 
parameters  may  be  defined  for  the  circular-to-coaxial  waveguide  junction. 
Figure  5b  illustrates  the  admittance  parameter  network  described  by 
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The  currents  in  (5.7)  are  related  to  the  incident  and  reflected  currents 
in  (5.2)  by  the  relations 

i.  -  it  ♦  i;  ;  v  =  it  -  r  (5.8 
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are  the  characteristic  admittances  for  the  equivalent  transmission  line 
representations  of  the  circular  and  coaxial  transmission  lines  respectively. 

Using  (5.7)  through  (5.11)  with  (5.2)  we  obtain  the  admittance 
parameters  in  terms  of  the  current  scattering  parameters,  which  are  given 
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D=  (1  -Cn)(  1  -C^)  -C12C21 


(5.12) 
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(5.14) 


(5.15) 


(5.16) 


As  in  the  case  of  the  previous  current  scattering  parameters,  the 
admittance  parameters  in  (5.12)  through  (5.15)  describe  a  lossless 
reciprocal  junction. 

The  simplest  form  of  (5.1)  occurs  when  the  circular  waveguide  contains 
no  "equivalent"  sources  or  propagating  modes,  thus  leaving  the  coaxial 
TMqo  (TEM)  mode  as  the  only  incident  mode  upon  the  circular-to-coaxial 
junction.  This  situation  corresponds  to  the  one-port  network  shown  in 
Figure  5c,  consisting  of  a  transmission  line  with  characteristic 


admittance,  Y_2»  given  in  (5.11)  terminated  by  the  capacitance, 

C  .  if1  4  Cz2  \  Y 
0  C2 


(5.17) 


where  oj  is  the  radian  frequency  and  C22  as  given  in  (5.6).  In  this 
one-port  case,  C22  is  the  reflection  coefficient  of  the  TEM  mode  current 
at  the  center  conductor  truncation.  The  capacitance,  CQ,  arises  from 
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the  energy  stored  in  the  evanescent  fields  of  the  p  =  a  and  p  =  b 
circular  waveguides  and  the  higher  order  fields  in  a  coaxial  waveguide. 

The  present  equivalent  capacitance  (which  is  for  a  hollow  center  conductor) 
is  compared  in  Figure  6  to  that  determined  through  a  variational  technique 
(which  considered  the  center  conductor  as  being  solid)  by  Risley  [13]  for 
a  British  3/4"  coaxial  line.  Due  to  the  difference  in  the  center  con¬ 
ductor  specification,  the  discrepancy  between  the  capacitances  is  not 
unexpected.  However,  our  value  of  the  equivalent  capacitance  reaches  a 
finite  limit  of  11.408  x  10  ^f  at  the  circular  TMqj  mode  cutoff  fre¬ 
quency.  It  cannot  be  determined  from  [13]  whether  the  variational  technique 
of  Risley  also  yields  a  finite  equivalent  capacitance  at  cutoff. 

The  electrical  effects  of  the  truncated  center  conductor  may  also 
be  expressed  in  terms  of  an  effective  lengthening  of  the  center  conductor 
as  shown  by  Marcuvitz  [12].  The  shift  of  the  open  circuit  reference  plane 
may  be  given  in  terms  of  the  present  parameters  as 


(5.18) 


Reference  plane  shifts  given  by  (5.18)  are  identical  to  those  given  by 
Marcuvitz  in  [12],  who  employed  the  same  transform  technique  as  used  in 
this  analysis. 

A  very  important  observation  to  make  at  this  point  is  that  the 
equivalent  one-port  network  for  the  circular-to-coaxial  waveguide  transition 
in  Figure  5C  may  be  simulated  by  terminating  the  "1"  port  of  the  two-port 
equivalent  circuit  in  Figure  5b  into  the  characteristic  admittance,  Y^, 
in  (5.10),  and  looking  into  the  "2"  port.  The  input  impedance  of  both 
systems  must  be  the  same, 


Thus  we  have  that 


Figure  6.  Equivalent  capacitance  of  a  coaxial  line  terminated  in  a  circular  waveguide. 

Solid  center  conductor,  Risley  [13],  hollow  center  conductor,  Equa.  (5.17). 
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which  is  valid  when  there  is  no  propagating  modes  in  the  circular  wave¬ 
guide  region  and  when  there  are  no  obstacles  too  near  the  junction  which 
would  cause  additional  reflections. 


6.  Electromagnetic  Penetration  into  the  end  of  a  Finite-Length  Thin  Cylinder 


Illustrated  in  Figure  1  is  a  finite-length  thin  (kb  «  1)  cylinder 
illuminated  by  a  uniform  plane  wave  incident  at  an  angle,  0^,  with  the 
z-axis.  Inside  and  recessed  from  the  end  of  this  hollow  outer  cylinder 
is  a  hollow  inner  cylinder.  We  assume  these  cylinders  to  have  infinitely 
thin  walls.  The  one  end  of  the  inner  cylinder  is  terminated  in  an 
arbitrary  load  admittance,  Y^.  Together  the  two  cylinder  comprise  a 
circular  waveguide  in  the  region,  -h  <z  <-h  and  a  coaxial  waveguide 

in  the  region,  -h  + <  z  <  -h  + 1  The  penetrating  current  on  the 

inner  wall  of  a  finite  length  hollow  thin  cylinder  illuminated  by  a 
uniform  plane  wave  was  derived  by  Chang,  et  al  [11],  and  may  be  written 
in  the  form  of  (3.1)  with 

<>b>0N  ’  -'V  ‘VOW  *  Cs(»-9i)I„(2h)TNW]E'  (6.1) 

The  reader  is  directed  to  the  above  reference  for  the  exact  definitions 
of  the  terms  in  (6.1).  Using  the  incident  current  modal  coefficients  in 
(6.1)  with  the  scattering  parameter  characterization  of  the  circular-to- 
coaxial  waveguide  transition  the  current  in  the  load  may  be  determined 
with  certain  restrictions.  First  the  analysis  in  [11]  requires  kb  «  1. 
Thus  all  of  the  circular  and  coaxial  modes  are  below  cutoff  except  of 
course  the  TMqq  (or  TEM)  coaxial  mode.  Secondly,  we  do  not  wish  to 
do  a  detailed  analysis  of  the  physical  make-up  of  the  load  and  its  effect 
on  the  reflections  and  mode  conversions  at  the  load,  thus  we  restrict 
the  analysis  by  assuming 

*c  >  b  (6.2) 
so  that  the  higher  order  coaxial  inodes  at  the  load  are  negligible  in 
comparison  to  the  TEM  mode.  And  by  requiring  that. 
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l  >  b  (6.3) 

w  " 

we  need  only  to  consider  the  TM^  mode  currents  in  the  circular  wave¬ 
guide  region,  the  effects  of  the  higher  order  ™on  mode  currents  being 
negligible  compared  to  the  effects  of  the  TM^  current.  Also  the 
restriction  in  (6.3)  permits  the  multiply  reflected  TM  mode  currents 
within  the  circular  waveguide  section  to  be  neglected  in  almost  all  cases. 
This  subject  is  addressed  mere  fully  later. 

The  restrictions  of  (6.2)  and  (6.3)  permit  the  use  of  the  equivalent 
two-port  networks  (C  or  Y  parameters  from  Section  5)  to  simulate  the 
circular-to-coaxial  waveguide  transition.  We  choose  the  admittance  para¬ 
meter  representation  described  in  (5.7)  over  the  current  scattering  para¬ 
meters  in  (5.2)  in  order  to  achieve  a  conventional  equivalent  circuit  in 
terms  of  currents  and  admittances  for  the  overall  system. 

Since  the  circular-to-coaxial  waveguide  transition  is  reciprocal 

(implying  Y ^  51  Y2l^’  we  maY  rcPrescut  this  junction  in  terms  of  a  simple 

tt  network.  The  overall  system  is  then  expressable  in  terms  of  the  circuit 

configuration  in  Figure  7a.  The  advantage  of  this  type  of  representation 

is  that  we  may  easily  introduce  equivalent  circuits  for  such  items  as; 

short  circuiting  plungers,  abrupt  changes  in  the  inner  or  outer  conductor 

radius,  multiple  loads,  etc.,  into  the  coaxial  line  section.  Note  also 

in  Figure  7a,  that  we  have  included  the  source  admittance  Yg  .  This 

"end"  admittance  term  corresponds  to  the  equivalent  admittance  seen  by  a 

TMqj  mode  current  incident  upon  the  open  end  of  the  cylinder.  From 

physical  considerations,  we  would  expect  Yg  to  contain  a  small  positive 

conductance  corresponding  to  the  effects  of  radiation  from  the  finite 

length  cylinder  and  for  | Y  j  to  be  quite  small.  In  this  section 

we  will  not  take  the  effects  of  Y  into  account  since  under  the  condition 

e 


u» 


Figure  7.  Equivalent  circuits  for  the  simulation  of  the  electromagnetic 
penetration  into  a  finite  length  coaxial  cylinder. 

(a)  Overall  system,  Note  in  general  the  coaxial  line  section  of 
length  £c,  may  include  additional  elements.  End  admittance, 

Y  ,  is  shown  here  but  neglected  in  analysis. 

(b)  Norton  equivalent  circuit. 
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in  (6.3),  the  input  admittance  to  the  circular  waveguide  at  the  circular- 
to-coaxial  waveguide  transition  is  essentially  Y^  ,  the  characteristic 
admittance  of  the  circular  waveguides  transmission  line  equivalent  given  in 
(5.10).  Essentially  we  are  neglecting  the  multiply  reflected  TM 
currents  within  the  circular  waveguide  section  by  assuming  the  circular 
waveguide  to  be  infinite  in  extent.  With  this  approximation  and  bearing 

g 

in  mind  that  (1^  )  in  (6.1)  is  an  incident  (rather  than  a  total) 
current  we  may  lump  the  entire  circuit  (including  the  n  network)  to  the 
left  of  the  coaxial  line  section  into  a  Norton  equivalent  circuit  as 
shown  in  Figure  7b.  The  equivalent  (total)  current  source  is  given  by: 


I 


N 


1  -2  Y  \  -y.  .  I 

_ 12  l  bl  w 

Y  ♦  Y  I  e 

*11  Yci/ 


(Ib  ')01 


(6.4) 


The  Norton  shunt  admittance  as  shown  in  Figure  7b  is  found  to  be 


YN  =  Y 


22 


Y  Y 
12  21 

Y  ♦  Y 
11  Cl 


(6.5) 


In  obtaining  (6.5)  we  have  again  neglected  the  multiple  reflections 
within  the  circular  waveguide  section  by  essentially  taking  the  circular 
waveguide  to  be  infinitely  long.  Comparing  (6.5)  with  (5.19)  we  find. 


Y 


N 


-i  at  C 


0 


(6.6) 


This  is  not  an  unexpected  result,  since  under  our  present  approximations, 
the  Norton  equivalent  shunt  admittance  has  the  same  definition  as  the 
discontinuity  admittance  (from  C^  )  as  used  in  Section  5. 

The  total  power  dissipated  within  the  coaxial  line  section  may  be 
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easily  calculated  from: 


(6.7)  does  not  require  any  intricate  circuit  analysis,  but  only  a  know- 

N  N 

ledge  of  the  Norton  equivalent  circuit  parameters,  I  in  (6.4)  and  Y 
in  (6.6),  and  the  input  admittance,  Y*n  ,  of  the  coaxial  line  section. 
The  maximum  power  dissipated  within  the  coaxial  line  section  would  occur 
when  a  conjugate  match  is  obtained  between  the  source  admittance  and  the 
coaxial  line  section  input  admittance,  i.e.,  Y*n  =  (YN)*  .  Linder  this 
conjugate  match  the  power  would  be  given  by, 


P 

m 


(6.8) 


where  G?  is  the  real  part  of  Y,,  .  We  immediately  observe  that  (6.8) 

predicts  an  infinitely  increasing  power  dissipation  for  a  conjugately 

matched  load  as  G*11  0  .  This  is  a  direct  consequence  of  neglecting  the 

end  admittance,  Y  ,  in  the  calculation  of  the  Norton  shunt  admittance, 
e 

N  N 

Y  ,  in  (6.5)  with  the  ultimate  result  of  Y  being  a  pure  imaginary 

N 

number.  Actually  had  the  exact  form  of  Y  been  used  in  (6.7)  the  power 
dissipated  in  the  coaxial  line  section  under  conjugately  matched  conditions 
would  have  been, 


P 

m 


(6.9) 


where  GN  is  the  real  part  of  Y^  .  The  maximum  possible  power  occurs 
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1 

in  N  N 

for  =  G  .  However  since  G  is  really  quite  small  we  may  confi¬ 

dently  use  (6.7)  and  (6.8)  to  calculate  the  power  in  the  coaxial  line 
section  except  for  conjugately  matched  cases  in  which  G*n  is  also  quite 
small. 

Utilizing  the  Smith  chart  as  a  means  of  illustrating  all  the  possible 

values  of  input  admittances  (normalized  to  Y  )  of  the  coaxial  line 

section,  it  is  obvious  that  conjugate  matching  only  occurs  along  the 

constant  susceptance  contour  Im [V^0]  =  (YN)*  .  And  it  can  be  shown 

from  (6.7)  that  the  contours  of  constant  power  are  circles  tangent  to 

the  point  Y^n  =  (Y^)*  .  In  both  cases  Y^  is  taken  to  be  as  given 

in  (6.6).  To  demonstrate  this  with  a  numerical  example,  we  consider  the 

situation  of  a  thin  cylinder  of  length  2h  =  1.30  m,  which  is  illuminated 

by  a  unit  (E*  =  1.0  V/m)  uniform  plane  wave  incident  at  an  angle 

0^  =  140°  with  the  z-axis  (see  Figure  1).  This  particular  incident  angle 

has  been  shown  in  [11]  to  provide  the  maximum  penetrating  current  into  the 

open  end  of  the  cylinder.  The  cylinder  investigated  here  has  a  finite 

outer  wall  thickness:  outer  radius  b  =  5.08  cm  and  inner  radius 

o 

b^  =  4.1275  cm  and  we  assume  the  penetrating  current  into  the  circular 
waveguide  region  to  be  equivalent  to  the  infinitely  thin  wall  case  where 
b  =  5.08  cm.  The  center  conductor  which  forms  a  coaxial  line  with  the 
inner  radius  of  the  larger  cylinder  is  recessed  by  the  distance, 
l  =  4.6  cm,  from  the  open  end  of  the  larger  cylinder  and  has  a  radius 
a  =  0.71374  cm.  In  Figure  8  are  shown  contours  of  constant  power  calculated 
from  (6.7)  as  a  function  of  the  normalized  input  admittance,  Y^/Y^  , 
of  the  coaxial  line  section.  Note,  in  Figure  8  we  have  expressed  the 
admittances  in  the  form  corresponding  to  the  conventional  time  convention, 
elut  .  Thus  the  coaxial  line  section  input  admittances  which  are  capacitive 


in  nature  are  in  the  upper  half  of  the  chart.  And  input  admittances  which 
are  inductive  are  in  the  lower  half  of  the  chart.  Along  the  constant 


susccptance  (normalized)  curve  B  =  -u)Cq/Y^ ,  =  -0.068  in  Figure  8,  the  power 
dissipated  in  the  coaxial  line  is  seen  to  increase  for  decreasing  values  of 
conductance.  In  Figure  9,  the  power  dissipated  in  the  conjugately  matched 
load  of  Figure  8  as  calculated  from  (6.8)  is  shown  as  a  function  of  the 
input  conductance  of  the  coaxial  line  section.  The  power  increases  indefinitely 
as  the  coaxial  line  input  conductance  is  made  smaller  and  smaller.  As 
previously  mentioned,  this  is  the  direct  result  of  assuming  the  circular 
wavguide  to  be  infinitely  long  for  the  determination  of  the  equivalent  Norton 
shunt  admittance.  In  realistic  terms,  however,  the  circular  waveguide  section 
is  finite  in  length  and  terminated  by  the  open  end  of  the  cylinder.  Radiation 
from  the  open  end  would  produce  a  real  component  in  the  equivalent  Norton 
shunt  admittance  and  thus  limit  the  maximum  power  that  could  be  dissipated 
within  the  coaxial  line  section.  The  maximum  possible  power  would  then  be 
given  by 


P  = 
m 


,N,2 


8G 


(6.10) 


N 

where  G  is  the  real  part  of  equivalent  Norton  shunt  admittance.  Further 
work  will  show  this  maximum  possible  power  to  be  independent  of  the  length 
of  the  circular  waveguide  (i.e.,  the  recession  of  the  inner  conductor). 
Although  for  deeper  recessions  of  the  inner  conductor,  conjugately  matched 
conditions  become  increasingly  difficult  to  obtain. 
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7.  Conclusions 

Through  the  independent  characterizations  of  the  various  elements 
comprising  the  finite  coaxial  line  (with  recessed  inner  conductor  and  load) 
illuminated  l>y  a  uniform  plane  wave  (figure  1),  we  have  constructed  a  simple 
equivalent  circuit  representation  to  yield  the  total  power  dissipated  within 
the  coaxial  line  section.  Under  the  mild  restriction  that  the  inner 
conductor  he  recessed  at  least  one  outer  cylinder  radius,  the  power  dissipated 
within  the  coaxial  line  section  under  most  situations  will  be  quite  small 
compared  to  the  maximum  possible  power  dissipation.  Higher  power  levels 
become  increasingly  difficult  to  obtain  due  to  the  more  exacting  specifica¬ 
tions  for  a  conjugately  matched  situation.  The  maximum  possible  power  is 
limited  by  the  real  part  of  the  equivalent  Norton  shunt  admittance,  which  is 
t lie  admittance  as  viewed  from  the  coaxial  line  side  of  the  circular-to-coaxial 
transition  circular  waveguide  section  terminated  by  the  open  end. 
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Appendix  A.  l;actori zation  of  M(ot)  and  N(a) 


Consider  the  kernel  in  (2.14), 


Vt!,» 


where  r,  is  given  in  (2.4).  Wc  wish  to  factorize  M {a)  into  the  form. 


N(a)  =  N  (a)N  (a) 


(A.  2) 


in  which  N+(a)  is  analytic  and  free  of  zeroes  in  the  upper  half  plane, 
-lin(k)  <  lm(a)  <  <*>  ,  and  N  (a)  is  analytic  and  free  of  zeroes  in  the 
lower  half  plane,  <  Im(a)  <  !m(k),  (sec  figure  4).  Actually  N (u ) 
in  (A. 1 )  is  characterized  only  by  simple  zeroes  at  u  =  ±  iy^,  defined 
in  (3. IK)  and  simple  poles  at  «  =  ±  iy^  ,  defined  in  (3.2).  There  are 
no  branch  cuts  or  other  irregularities.  Since  both  .1  ^ (r,  )  and  d^K^) 
are  analytic  functions  of  cx,  they  may  each  be  expressed  in  terms  of 
infinite  products  of  their  respective  zeroes.  following  Whittaker  and 
Watson  1 1  (» | , (A. 1 )  may  be  expressed  in  terms  of  the  infinite  product, 


-  *  f-H1 

.  pj-  1,  V - laaii. 


^bn) 


(A.  3) 


where, 


J0(ka) 


0  1/  ^()(kb) 


(A. 4) 
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After  the  introduction  of  a  couple  of  exponential  factors  in  the  infinite 
product  (which  provide  convergent  individual  products  but  cancel  in  (A. 3)), 
we  may  factorize  (A. 3)  into  the  form  of  (A. 2)  where. 


N+(cxJ  =  Nq  n 
n=l 


(A. 5) 


and  N  (a)  has  a  similar  form  except  that  (a)  is  replaced  by  (-a)  in 
all  terms. 

To  speed  the  convergence  of  (A. 5)  we  follow  the  procedure  of 
Olson  [17],  which  utilizes  the  relationship  from  Mittra  and  Lee  [14], 


II  e 
n=  1 


(A. 6) 


where  T(x)  is  the  gamma  function  and  y  =  0.5772156...  .  Choosing  the 

appropriate  values  of  A  and  B  in  each  product  term  of  (A. 5)  and  (A. 6)  and 
incorporating  (A. 7)  with  these  respective  values  in  (A. 5)  and  (A. 6)  results 
in  very  quickly  convergent  infinite  product  expressions  for  N+(a)  and 
N  (a).  For  the  former  we  have. 


N+(a)  =  N0exp[-iy — (1-v)] 


(A- 7) 


where  v  =  a/b  and  y^^  and  ^bn^  are  t*ie  asymPtot*c  forms  as  n  -*■  00 
of  Yan  in  (3.18)  and  y^  in  (3.2),  respectively,  and  are  given  by. 


,(°)  , 

'  an 


Y„V  =  ("  -  1/4)  £ 


bn 


=  (n  -  1/4) 


(A. 8) 
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Next  consider  the  kernel,  M(a),  in  (2.14), 

M(a)  =  n|.J()(ca)Y()(r,b)  -  .]()(.\h)Y()(r,a )  |  (A.  9) 

We  wish  to  factor  ice  N(u)  into  the  form. 


M(ct)  =  M+  (aJM  (a) 


(A.IOj 


in  which  M  (n)  and  M_(tx)  arc  analytic  and  free  of  zeroes  in  the  upper 
and  lower  half  planes,  respectively.  Since  M(u)  in  (A. 9)  is  an 
analytic  function  of  u,  possessing  only  simple  zeroes  at  u  =  i  iy ,  , 
defined  in  (3.X),  we  may  follow  (he  same  steps  that  led  to  (A. 11)  for 
N+(ix).  And  we  get, 


ctb 


l+(a)  =  M()cxp{  iy  —  (1-v) } 


II 

rn=  1 


where, 


M()  =  •/^CJ(ka)Yo(kh)  -  J0(kh)Y()(ka) 

and , 


(A. 12) 


(Oj  _  Jim 
fcm  “  (b-aj 


(A.  13) 


y^  is  the  asymptotic  form  of  y  in  (3.H)  as  n  -*■  00  .  We  also  find 
cm  cm 

that  M  (u)  is  of  the  same  form  as  M+(u)  in  (A.l<>)  except  (a)  is  replaced 
by  (-a) 


1 1 J 


r>  o 


The  introduction  of  complementary  entire  functions  into  N+(u)  and 
N  (a)  (i.e.,  these  functions  cancel  in  the  product  N(a)  =  N  (<xjN  (u) ) 

docs  not  alter  the  regions  of  analyticity  of  the  respective  factors.  The 
same  applies  to  the  factors  M  (a)  and  M  (a)  .  Thus  in  order  to  obtain 
the  desired  algebraic  asymptotic  behavior  of  the  product  M+ft)N+(  f|  as 
|u|  '  *■  in  the  upper  half  plane,  we  may  introduce  the  entire  functions 
expt-i ~  l-y(l-v)  +  v  tn  (v) | }  and  exp{-i  ~  ( 1-v) |y  +  ^n(l-v)]}  into 
N+ft)  and  M  (u),  respectively.  The  reciprocals  of  these  functions  are 
introduced  into  the  respective  minus  functions  also,  so  that  the  original 
kernels,  N(u)  and  M(u),  are  still  obtained  from  the  product  of  their 
respective  plus  and  minus  factors.  With  the  above  modification,  we  get 
the  final  forms  of  the  plus  factors, 


<xl) 

N  ('i)  =  N  .exp |  -  i  -  v  c n  (  v )  I 

4-  II  ti 


/  . 

I-' 
'.  1 


-  i 


n  /  V  an  / 


1  + 


(0) 

bn 


('  ■  ""  l1  *  i1  * 

an 


(A. 14) 


and 


1  ♦ 


M+(u) 


M()exp(-i  ~  ( l-v)Hn(l-v)}  — 


i  Y 

cm  i 


(l  -  i^p(l-  v))  m=  1  / 1 


1Y 


(0) 

cm  , 


(A. IS) 


The  corresponding  minus  factors  arc  similar  in  form  to  the  plus  factors 
except  fx)  is  replaced  by  (-a)  throughout,  i.e.. 


5  1 


N_  («)  =  N+{-(i)  (A.  10) 

and 

M  (a)  =  M  (-0  (A.  17) 

The  desired  algebraic  asymptotic  behavior  of  the  product  M  (a)N  (uj 
is  given  by, 

M+(u)N+tu)  'v>  a  ;  as  |'x(  ►  «>  in  the  upper  half  plane  ( A .  1 S ) 

It  is  worthwhile  to  mention  that  N +(u)  is  an  analytic  function  of 
ot  even  in  the  lower  half  plane  except  at  the  poles  located  at  a  =  -iy^, 

n  =  1,2,5 .  And  [N+(a)]  1  is  also  an  analytic  function  of  a  in  the 

lower  half  plane  except  at  the  poles  located  at  u  =  -iy  ^  =  1,2,5,... 

Thus  evaluating  N+ (a)  or  [N+(a)J  1  in  the  lower  half  plane  poses  no  real 
difficulties.  liven  evaluating  a  plus  factor  at  one  of  t lie  fore-mentioned 
poles  does  not  yield  an  unphvsical  result  in  the  solutions  for  the  current. 
This  situation  is  discussed  in  Appendix  B.  Similar  arguments  apply  to  the 
evaluation  of  M+(<x)  and  [M  (ct )  j  ^  in  the  lower  half  plane. 

finally,  when  using  (A.  14)  to  obtain  N+(ci)  or  (A.  15)  to  find  M+  (a) , 
it  is  necessary  to  compute  (depending  on  the  magnitude  of  a)  only  a  limited 
number  of  products.  This  is  due  to  the  rapidly  convergent  nature  of  the 
infinite  products  in  both  (A. 14)  and  (A. 15). 


Appendix  If . _ I  tie  "Wei  1  -dimpled"  Modes 


In  the  ea  I  cu  1  ;i  t i on  of  the  residue  terms  for  the  integrals  in  (4.19) 

resulting  in  the  expressions  for  I^(z);  (4.20)  for  z-_0  and  (4.21  J  for 

,  wo  have  implicitly  assumed  (and  correctly  so)  that  the  poles  of  the 

integrands  were  not  coincident,  i.e.,  they  were  all  separate  first  order 

poles.  this  leads  to  the  observation  of  the  questionable  terms: 

(  y .  -  y  ,.)  1  in  (1.2U)  and  (y  -  y,,.)  1  in  (4.21),  which  seemingly 

bn  cM  cm  b.Ni 

make  the  current  "blow  up"  if  n  =  p  and  p  =  p,.  ,  respectively. 

n  cM  cm  N 

fortunately,  this  does  not  happen  since  the  terms  M  (-iy,  )  and  M  ( -  i  y .  ,, ) 

11  +  bn  +  bN 

are  zero  at  p  =  p  ,,  and  p  -  m,  .  Actually  both  of  the  above  situa- 
n  cM  cm  N 

tions  will  occur  at  the  same  time,  that  is,  whenever  a  zero  (p  or  p,  ) 

n  N 

of  d  (x)  coincides  with  a  zero  (p  or  p^ )  of  |.J()(vx)Y()(x) -J(|(x)  Y()(vx)  J 

We  denote  these  coincident  zeroes  as  p  =  p  ..  and  p,,  =  P  „  ,  which 

q  cT  Q  cp 

differ  only  in  their  individual  relationships  to  the  incident  (capital 
letters  in  subscript,  I’  and  Q)  and  scattered  (lower  case  letters  in  sub¬ 
script,  p  and  q)  modes. 

The  expressions  for  l^(z)  in  Section  4,  ( 4 . 2 1* )  for  z^O  and  (4.21) 
for  Zj(>  ,  are  perfectly  valid  in  the  "wol 1 -coupled"  mode  case  if  the 
following  limits  arc  observed. 
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(If.l)  and  (If.  2)  contain  the  terms. 


a'M\p=iV, 


VVpc!’* 


W 


(yp„  ! 


N(a)  =  v 


a>‘Ycp=lYbQ 


J ,  Co  ) 

1  cj> 


VV<W 

W 


The  substitution  of  the  limits  in  f  H  -  1 J  through  (11.11  into  (4  .Jo)  and 

(4.21)  whenever  coincident  zeroes  occur  will  yield  the  correct  result  for 

l^(z)  in  the  "wol 1 -coupled"  mode  ease. 

Mathematically  these  "well-coupled”  modes  occur  whenever  an  eigenvalue 

(My,  or  My,,.)  of  the  circular  waveguide  corresponds  exnctlv  to  an 
bq  bQ 

eigenvalue  (My^  or  ±iy  )  of  the  coaxial  waveguide.  Physically  the 
"well-coupled"  modes  amount  to  a  circular  waveguide  mode  (cither  incident 
from  z<0  or  scattered  due  to  the  coaxial  source  for  zMl)  having  a  zero 
z-directcd  electric  field  at  p=a  (the  center  conductors  radial  position). 
The  particular  "wel 1 -coupl ed"  circular  and  coaxial  waveguide  modes  then 
match  fairly  well,  although  not  exactly,  at  the  transition  at  z=0  .  We 
note  that  there  is  nothing  discontinuous  in  the  "wel 1 -coupled"  mode  pro¬ 
cess  and  the  closely  matched  situation  between  the  circular  and  coaxial 
waveguide  modes  is  approached  from  both  sides  of  the  "wel 1 -coupl ed"  con¬ 
dition.  The  lowest  order  modes  for  which  this  effect  may  be  observed 
are  the  circular  I'M,.,  and  the  coaxial  TM  modes.  In  this  situation, 

i)  =p  =r».r>2l)l .  The  coaxial  TM,,,,  lor  TPM)  mode,  however,  will  never 

2  C  1  (HI 

he  "wel I -coupled"  to  any  of  the  circular  waveguide  modes. 


